Search for: All records

Based on a generalized local Kolmogorov–Hill equation expressing the evolution of kinetic energy integrated over spheres of size$\ell$in the inertial range of fluid turbulence, we examine a possible definition of entropy and entropy generation for turbulence. Its measurement from direct numerical simulations in isotropic turbulence leads to confirmation of the validity of the fluctuation relation (FR) from non-equilibrium thermodynamics in the inertial range of turbulent flows. Specifically, the ratio of probability densities of forward and inverse cascade at scale$\ell$is shown to follow exponential behaviour with the entropy generation rate if the latter is defined by including an appropriately defined notion of ‘temperature of turbulence’ proportional to the kinetic energy at scale$\ell$.

 

The interplay between viscoelasticity and inertia in dilute polymer solutions at high deformation rates can result in inertioelastic instabilities. The nonlinear evolution of these instabilities generates a state of turbulence with significantly different spatiotemporal features compared to Newtonian turbulence, termed elastoinertial turbulence (EIT). We ex- plore EIT by studying the dynamics of a submerged planar jet of a dilute aqueous polymer solution injected into a quiescent tank of water using a combination of schlieren imaging and laser Doppler velocimetry (LDV). We show how fluid elasticity has a nonmonotonic effect on the jet stability depending on its magnitude, creating two distinct regimes in which elastic effects can either destabilize or stabilize the jet. In agreement with linear stability analyses of viscoelastic jets, an inertioelastic shear-layer instability emerges near the edge of the jet for small levels of elasticity, independent of bulk undulations in the fluid column. The growth of this disturbance mode destabilizes the flow, resulting in a turbulence transition at lower Reynolds numbers and closer to the nozzle compared to the conditions required for the transition to turbulence in a Newtonian jet. Increasing the fluid elasticity merges the shear-layer instability into a bulk instability of the jet column. In this regime, elastic tensile stresses generated in the shear layer act as an “elastic membrane” that partially stabilizes the flow, retarding the transition to turbulence to higher levels of inertia and greater distances from the nozzle. In the fully turbulent state far from the nozzle, planar viscoelastic jets exhibit unique spatiotemporal features associated with EIT. The time-averaged angle of jet spreading, an Eulerian measure of the degree of entrainment, and the centerline velocity of the jets both evolve self-similarly with distance from the nozzle. The autocovariance of the schlieren images in the fully turbulent region of the jets shows coherent structures that are elongated in the streamwise direction, consistent with the suppression of streamwise vortices by elastic stresses. These coherent structures give a higher spectral energy to small frequency modes in EIT characterized by LDV measurements of the velocity fluctuations at the jet centerline. Finally, our LDV measurements reveal a frequency spectrum characterized by a −3 power-law exponent, different from the well-known −5/3 power-law exponent characteristic of Newtonian turbulence. 

Prior mathematical work of Constantin & Iyer ( Commun. Pure Appl. Maths , vol. 61, 2008, pp. 330–345; Ann. Appl. Probab. , vol. 21, 2011, pp. 1466–1492) has shown that incompressible Navier–Stokes solutions possess infinitely many stochastic Lagrangian conservation laws for vorticity, backward in time, which generalize the invariants of Cauchy ( Sciences mathématiques et physique , vol. I, 1815, pp. 33–73) for smooth Euler solutions. We reformulate this theory for the case of wall-bounded flows by appealing to the Kuz'min ( Phys. Lett. A , vol. 96, 1983, pp. 88–90)–Oseledets ( Russ. Math. Surv. , vol. 44, 1989, p. 210) representation of Navier–Stokes dynamics, in terms of the vortex-momentum density associated to a continuous distribution of infinitesimal vortex rings. The Constantin–Iyer theory provides an exact representation for vorticity at any interior point as an average over stochastic vorticity contributions transported from the wall. We point out relations of this Lagrangian formulation with the Eulerian theory of Lighthill (Boundary layer theory. In Laminar Boundary Layers (ed. L. Rosenhead), 1963, pp. 46–113)–Morton ( Geophys. Astrophys. Fluid Dyn. , vol. 28, 1984, pp. 277–308) for vorticity generation at solid walls, and also with a statistical result of Taylor ( Proc. R. Soc. Lond. A , vol. 135, 1932, pp. 685–702)–Huggins ( J. Low Temp. Phys. , vol. 96, 1994, pp. 317–346), which connects dissipative drag with organized cross-stream motion of vorticity and which is closely analogous to the ‘Josephson–Anderson relation’ for quantum superfluids. We elaborate a Monte Carlo numerical Lagrangian scheme to calculate the stochastic Cauchy invariants and their statistics, given the Eulerian space–time velocity field. The method is validated using an online database of a turbulent channel-flow simulation (Graham et al. , J. Turbul. , vol. 17, 2016, pp. 181–215), where conservation of the mean Cauchy invariant is verified for two selected buffer-layer events corresponding to an ‘ejection’ and a ‘sweep’. The variances of the stochastic Cauchy invariants grow exponentially backward in time, however, revealing Lagrangian chaos of the stochastic trajectories undergoing both fluid advection and viscous diffusion. 

We use an online database of a turbulent channel-flow simulation at $Re_\tau =1000$ (Graham et al. J. Turbul. , vol. 17, issue 2, 2016, pp. 181–215) to determine the origin of vorticity in the near-wall buffer layer. Following an experimental study of Sheng et al. ( J. Fluid Mech. , vol. 633, 2009, pp.17–60), we identify typical ‘ejection’ and ‘sweep’ events in the buffer layer by local minima/maxima of the wall stress. In contrast to their conjecture, however, we find that vortex lifting from the wall is not a discrete event requiring $\sim$ 1 viscous time and $\sim$ 10 wall units, but is instead a distributed process over a space–time region at least $1\sim 2$ orders of magnitude larger in extent. To reach this conclusion, we exploit a rigorous mathematical theory of vorticity dynamics for Navier–Stokes solutions, in terms of stochastic Lagrangian flows and stochastic Cauchy invariants, conserved on average backward in time. This theory yields exact expressions for vorticity inside the flow domain in terms of vorticity at the wall, as transported by viscous diffusion and by nonlinear advection, stretching and rotation. We show that Lagrangian chaos observed in the buffer layer can be reconciled with saturated vorticity magnitude by ‘virtual reconnection’: although the Eulerian vorticity field in the viscous sublayer has a single sign of spanwise component, opposite signs of Lagrangian vorticity evolve by rotation and cancel by viscous destruction. Our analysis reveals many unifying features of classical fluids and quantum superfluids. We argue that ‘bundles’ of quantized vortices in superfluid turbulence will also exhibit stochastic Lagrangian dynamics and satisfy stochastic conservation laws resulting from particle relabelling symmetry. 

Transition from laminar to turbulent flow occurring over a smooth surface is a particularly important route to chaos in fluid dynamics. It often occurs via sporadic inception of spatially localized patches (spots) of turbulence that grow and merge downstream to become the fully turbulent boundary layer. A long-standing question has been whether these incipient spots already contain properties of high-Reynolds-number, developed turbulence. In this study, the question is posed for geometric scaling properties of the interface separating turbulence within the spots from the outer flow. For high-Reynolds-number turbulence, such interfaces are known to display fractal scaling laws with a dimension$D\approx 7/3$, where the 1/3 excess exponent above 2 (smooth surfaces) follows from Kolmogorov scaling of velocity fluctuations. The data used in this study are from a direct numerical simulation, and the spot boundaries (interfaces) are determined by using an unsupervised machine-learning method that can identify such interfaces without setting arbitrary thresholds. Wide separation between small and large scales during transition is provided by the large range of spot volumes, enabling accurate measurements of the volume–area fractal scaling exponent. Measurements show a dimension of$D=2.36\pm 0.03$over almost 5 decades of spot volume, i.e., trends fully consistent with high-Reynolds-number turbulence. Additional observations pertaining to the dependence on height above the surface are also presented. Results provide evidence that turbulent spots exhibit high-Reynolds-number fractal-scaling properties already during early transitional and nonisotropic stages of the flow evolution.

 

We consider the problem of formulating perturbative expansions of the conformation tensor, which is a positive definite tensor representing polymer deformation in viscoelastic flows. The classical approach does not explicitly take into account that the perturbed tensor must remain positive definite – a fact that has important physical implications, e.g. extensions and compressions are represented similarly to within a negative sign, when physically the former are unbounded and the latter are bounded from below. Mathematically, the classical approach assumes that the underlying geometry is Euclidean, and this assumption is not satisfied by the manifold of positive definite tensors. We provide an alternative formulation that retains the conveniences of classical perturbation methods used for generating linear and weakly nonlinear expansions, but also provides a clear physical interpretation and a mathematical basis for analysis. The approach is based on treating a perturbation as a sequence of successively smaller deformations of the polymer. Each deformation is modelled explicitly using geodesics on the manifold of positive definite tensors. Using geodesics, and associated geodesic distances, is the natural way to model perturbations to positive definite tensors because it is consistent with the manifold geometry. Approximations of the geodesics can then be used to reduce the total deformation to a series expansion in the small perturbation limit. We illustrate our approach using direct numerical simulations of the nonlinear evolution of Tollmien–Schlichting waves. 

This work introduces a mathematical approach to analysing the polymer dynamics in turbulent viscoelastic flows that uses a new geometric decomposition of the conformation tensor, along with associated scalar measures of the polymer fluctuations. The approach circumvents an inherent difficulty in traditional Reynolds decompositions of the conformation tensor: the fluctuating tensor fields are not positive definite and so do not retain the physical meaning of the tensor. The geometric decomposition of the conformation tensor yields both mean and fluctuating tensor fields that are positive definite. The fluctuating tensor in the present decomposition has a clear physical interpretation as a polymer deformation relative to the mean configuration. Scalar measures of this fluctuating conformation tensor are developed based on the non-Euclidean geometry of the set of positive definite tensors. Drag-reduced viscoelastic turbulent channel flow is then used an example case study. The conformation tensor field, obtained using direct numerical simulations, is analysed using the proposed framework.